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The  general  theory  of  stationary  spatial  models 
is  developed:  namely  MA,  moving  average;  AR,  auto¬ 
regressive;  and  ARMA,  autoregressive  moving  aver¬ 
age  processes.  As  compared  to  the  time  series  in 
m  dimensions,  spatial  models  may  be  one-sided, 
two-sided,  or  mixed.  Free  use  is  made  of  the 
previous  results  of  Aroian  and  his  associates  in 
time  series  in  m  dimensions.  The  main  theoreti¬ 
cal  properties  of  the  models  in  the  univariate 
case  are  established.  The  multivariate  case  is 
even  more  important  than  the  univariate.  Estima¬ 
tion  by  minimum  variance  and  simulation  of  the 
models  are  included. 
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RODUCTION 


The  results  of  time  series  in  ra  dimensions  by 
Aroian  and  his  coauthors  are  used  to  establish 
the  results  of  spatial  models  in  m  dimensions. 

If  m=l  the  results  apply  to  events  on  a  line  such 
as  a  river  at  a  particular  time;  for  ra=2,  the 
events  are  those  in  the  plane  such  as  ecological 
distribution  of  a  plant,  the  average  rainfall  for 
the  plane  after  a  storm  is  completed;  for  m=3, 
pollution  is  space  or  distribution  of  ^mineral  in 
a  mine. 

Important  assumptions  are  outlined:  the  charac¬ 
teristic  of  an  event  in  space  is  given  by 

Z  ,X=  (X_ ,X-,  . . .  ,x  ),  — 
x  l  z  m 

x— i— (x. — — Z  ,  . . • , x-fc) . 

I  I  2.  z  mm 

Weak  stationarity  is  assumed  in  space  as  a  mini¬ 
mum  assumption: 

2  2 

y  =e (z  ) =0,  a  -E (z  -y  ) 

2S  X  Z  X  Z 

E(ax)=o,  o2>o,  e*«K<V*-i-V2/0*' 

«.=  <«.,  ,JL . I  ) . 

i  2.  m 

All  second  order  moments  exist.  Note  x  may  be  in 
any  coordinate  system,  and  t  may  be  plus  or  minus; 
the  results  are  in  m  dimensions,  only  the  time 
coordinate  has  been  dropped  from  consideration. 
Although  time  is  a  variable,  it  is  not  spatial, 
so  new  theory  must  be  developed. 

Two  results,  in  general,  follow  from  time  ser¬ 
ies.  If  m=l ,  one-sided  spatial  models  are  cover¬ 
ed  by  Box  and  Jenkins  (1976)  if  the  variable  t  in 
their  models  is  replaced  by  X*  Isotropic  models 
in  space  m=2 ,  where  x  is  the  radius  of  a  circle, 
are  models  of  w*l,  time  series  in  m  variables, 
and  are  discussed  briefly  in  a  later  section. 

2.  MA,  AR,  AND  ARMA  MODELS 

The  two-sided  theoretical  spatial  MA  model  is 
defined: 

z  .  ?  Uj  a  ,  -«kz<»,  i/f  *1, 

x  n—  n  x-n  0  (2.1) 

n  -  (n  ,  n, . . . ,n  ) ,  £  •  £  ...  I  _ 

12  m  n*-«  n  ■-**  n. 

m  2  1 

a  is  an  i.i.d,  variable  with  y  "0,  a^> 0,  indepen¬ 
dent  of  z  unless  z  or  z  0  involves  a#  or 

X  X  X—  X  X  X—  Mt 


and  Ea  z  D*0,  unless  £»0.  /**' 

x  x-fc 

More  usually  n  is  finite: 

q 

z  *  L  a  ,  *  «1  (2.2) 

x  n»-*°  n  x+n  0 

an  MA  model  of  spatial  order  p.+q.  i"  «*ch  spat¬ 
ial  variable  x^  l<i<m.  If  -p^n^the  spatial 

model  is  two-sided  in  m;  if  ~p^£n£-q^  or  O^nSq^ 

for  all  i  it  is  one-sided  in  m.  A  model  may  be 
two-sided  for  certain  x.  and  one-sided  for  other 
x. ;  such  a  case  is  called  a  mixed  model. 

As  examples:  for  m*l, 

v*i  “x-i^-i  vi^'xA'^x  |  (2-3) 

m=l ,  two-sided  of  order  two  in  each  variable.  If 
^  2*0,  the  model  is  one-sided.  If 

^2^0,  ^  2^0,  it  is  mixed. 


Z*l,x2*^01  ax^,x2-l+^10  ^j^-l ,x2+^0-lax^,x2+l 

♦*-10  \1+l.x2+ax1(x2  (2-4) 

of  spatial  order  two  in  each  variable,  a  two-sided 
model;  one-sided  model  if  ip -»ip  , 0*°f  and  mixed 

vr0'  *-ifc/0- 

The  two-sided  theoretical  spatial  AR  model  is 
defined: 

z  «=  I  $  z  +a  ,0  «0,  -«kz«»  .  (2.5) 

x  n— 00  n  x+n  x  0 

Usually  n  is  finite  -p<n<q,  of  spatial  order  p> 
q.  in  each  variable  ,  l£i£m.  It  may  be  two- 

sided,  one-sided  or  mixed  as  in  the  MA  model. 
Zx“Vx-l+*-lZx+l+*2Zx-2+*-2Zx+2+ax  (2*6) 

?X1  »X2+^^Zxl  ,x2+* 

+$  +a  »  (2.7) 

-10  x1+l,x2 

The  theoretical  ARMA  model  for  time  series  in 
m  dimensions  is,  Voss  et  al  (1980) : 

q  r  vs 

Zx,t*  n*-p  k-1  ^n,k  Zx+n,t-k  n--u  k*l  ^n,k 

ax+n,t-k+ax,t’  *00“°'  600-°  (2-8) 

(r,s)  in  the  temporal  domain,  q+p  and  u+v  in  each 
spatial  variable.  The  general  case  would  be 
-«<n<«»,  -“Xt<®,  The  corresponding  two-sided  ARMA 
spatial  model  is 

z  -  ?  <p  z  -  I  9  a  ,+a  ,  <f>  “0  -0.  (2.9) 

x  n*-p  Tn  x+n  n«-u  n  x+n  x  0  0 

Examples  for  m-1  and  2  respectively  are: 

Zx"*l  ZX-1**2  Zx-2+-l  Zx+l+*-2Zx+2~*  1  *x-l 

-eiax-re2ax-2-e-lax+r®-2ax+2+ax  <2*10) 

ARMA  model  two-sided  p-q*2 ,  m*l . 


*x1,x2'v0l‘x1,x2-rv10*xx-l,x2^0-l-x1,x2+l 

4^-10Zx1+l,x2"®01ax1  ,x2-l  ®10axx-l,x2 


80 
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’00— r'xj  ,x  ,*r°-I<'‘'x1*l,x  ,*',x|  .x  ,  •  (2-11) 

AKMA  model,  two-s i di‘U  m  . 

The  corresponding  two-sided  MA  (AR)  model  is 
found  by  setting  $+1*<fr+2«$+  3«=0  ( 0+1»O+2*0^3«O) 

and  the  one-sided  ARMA  model  by  setting 

<t> .  =0 ,  or  A  .-0  and  0.*O  or  0  .*0. 
i  -l  x  -i 

There  may  be  a  mixture  such  as  an  AFMA  model  with 
one-sided  AR  or  MA  while  the  other  part  is  two- 
sided. 

3.  MA  MODELS 
The  simplest  MA  model 

z  =-6,a  -  9  .a  ,  +  a  (3.1) 

x  1  x— 1  -1  x+1  x 

is  obtained  from  (2.3)  by  replacing  ^'s  by  -0's 
and  setting  l|>_2=-0_2,  ^2=-e2,  and  ®2=0_23S°#  n-1. 

For  this  model 

°H2(i+ei+0-i>  (3-2> 

P1=-(e1+e-1)/(i+02+ef;i) , 

P2"(eie-l,/a+0X+6-l)'  P1=P-1'  p2=P-2' 
and  all  other  p£=0,  i>2. 

If  (3.1)  reduces  to  the  one-sided 

model  in  t;  and  also  if  6^=0,  0  ^0.  Given 
(01#  ®2^»  are  g*ven  ky  (3.2).  Conversely 

given  a  sample  estimate  of  (p^p^,  (§  ,0  ) 

may  be  found  from  the  set  replacing  population 
values  by  estimates: 

2  2 

1 ,0^.0*  _  t  r\  .  r\  »  /  -  nr* 


1+ei+6-l=  ‘t8l+6-l)/Pi 


616-1/P2!  and 


Pj  =  -(1/01+1/0_1)P2,  (3.3) 

Set  0^=u+v,  0  1=u-v,  then 

1+2(u2+v2)«  -2u/p^  =(u2-vVp2/  (3.4) 

whose  solution  involves  the  intersection  of  a 
circle  and  hyperbola.  This  may  lead  to  four 
possible  sets,  but  the  condition 

I  ei  I  +  l  9-i  I  <  1  <3*5) 

limits  the  results  to  one  set. 

If  0  =6  _,  p-0,  and  if  0,  or  6  *0,  p  «0. 

1-11  1  -12 

Table  1  lists  the  values  of  {p^rP^}  given 
and  conversely.  The  manipulations  of  (3.4)  show 
|P2I  <  1/3,  and  (l+2p2)2  2  4p2. 

Only  the  values  of  {p^Pj}  for  -1  £  6  $  1,  and 

0^2  <  ®  are  tabled,  since  the  remaining  values  of 
{02'®2^  f°un<i  from  the  skew  symmetry 

implied  by  (3.3) . 

The  characteristic  equation  of  (3.1)  is 
1-02  Bx“e_2  BJ(1f  and  the  corresponding  AR  repre¬ 
sentation  of  the  MA  model  is: 

ax  -  dlo  <ei  V8-lBx1,d  V  WVx-l 

+  5lzx+l  +  81Zx~2  +  28l8-l*x  +  8-1*x-2  +  ”• 


or  a  «  *nz  ♦»,*  .4*  *  ,  +  ...  (3.6) 

x  Oxl  x-J  -1  x+1 

Not  o  that  n,.  j  1.  arid  It  «*xt»dri  and  the  am  Ian 
0  n 

converges  for  the  values  ot  1 0^  |  +  1 0  |  ‘  1. 

Given  a  sample  of  n,  {0^,0  ^ }  is  estimated 
from  (3.3)  using  sample  estimates  (p^,P2)  and  the 

methods  of  solution  already  indicated.  The 
~  A  2  *2 

approximate  variance  of  0  is  og  -  (1-0  )/n#  and 

confidence  intervals  may  be  obtained  if  it  is 
assumed  that  the  errors  are  distributed  normally. 
Another  method  of  estimation  is  to  vary  0,  and 
choose  the  6  which  minimizes  the  variance  of  the 
error  of  prediction,  e-z^-z^. 

Another  simple  MA  model  m»2 ,  of  the  first 
order  is: 

Zx^,x2  ^01  ax^ ,x2-l"*^10ax^-l,X2  'x2+*- 

“e-10ax1+l,x2+ax1,x2  (3.7) 

The  characteristic  function  is: 

1-801  Bx2-®10  Bx1-80-1  B;2-8-10  B;;  '  (3-8) 

the  values  of  a2,  Pox.P1q.Po2*‘,20'p11*p-U  are: 

2  2  2  2  2  2 
°z  -  V1+0Ol*0O-l+0lO+0-lO>' 

P01  "  P0-l  ■  -(801+80-1>/Oz  \  ' 

P10  =  P-10  '  -«V8-10,/0z  CI2  ' 

P02  =  P0-2  “  t80180-ll/a2°a2 

P20  *  p-20  *  (810e-10,/tJz  °a 

pu  "  p-i-i  *  <W0o-i9-io,/oz 


P-ll“  Pl-1 


te„,  e  ,  +e„  ,e,jo2<r2 

01  -10  0-1  10  z  a 


1+8oi+8o-i+8L+0!io'  -(6oi+9o-i)poi  *  •(610+8-10> 
pIo  “  (6O101O+0O-18-1O,PU  -  (0O19-1O+0O-181O> 

P-ll  *  l0Ol0O-l,PO2  *  (01O0-1O)P2O  •  (3.10) 

The  corresponding  AR  model  is  given  by  the 
inversion  of  the  characteristic  function  in  the 
usual  way.  Note  |e01 k I e10  N ®0_1 I + I e_io I  <1. 
Given  a  set  ( pQ1  *P10»P11  *3^  the  proper  set 
o*  equations  from  (3.10)  are  used  to  estimate 
(801'®10  '®0-l'®-l  o* •  The  aPPtoximate  variances 

of  the  set  of  the  0*s  may  be  found  using  the 
methods  suggested  in  Aroian  and  Taneja  (1980) , 
Perry  and  Aroian  (1979) ,  and  Aroian  and  Schmee 
(1980) .  The  variance  of  the  error  of  predictions 
may  be  minimized  by  changing  the  vector  of 
estimates  6  until  a  minimum  is  found  for  this 
variance,  see  Aroian  and  Taneja  (1980). 

Simulation  of  MA(1,1)  Model 

I*t  0^-. 2,  0^— 72  in  (3.1) ,  a^'s  being  dis¬ 
tributed  as  N (0, 1) .  This  model  is  simulated  with 
100  observations.  First  random  shocks  *x'*  are 


generated,  and  z  *s  are  found  from  (3.1).  The 
estimated  correlation  coefficients  r. ,  r  are  ob¬ 
tained.  The  {00  is  found  from  T3.3;-(3.4) 

using  r^,  r^.  Minimum  error  prediction  of  0^  and 

6  ^  are  also  found  by  using  a^'s  obtained  from 

z^'s  by  minimizing  e^-tz-z^),  this  estimate  is 

{01,0  i}={ . 330, -.340} ,  note,  that  given  z^a^  * 

z  +0  a  +0  a  With  twenty-five  simulation 

x  1  x— 1  -1  x+1 

runs  it  is  found  that  $i=.284,  F  ^=-.239,  while 

minimum  variance  estimates  are  0  =.285,  9  =-.237. 

A  2  ^2  *2*2  1 

From  aG  Ml-9  )/n,  0A  =Oa  =.0384.  The  approximate 
o  0-  o-l 

1  a  a  -l 

formula  for  the  covariance  J (0^ , ©— l) n 

or  -.00784  for  the  simulated  case  versus  the 
actual  of  -.0032. 

4.  AR  MODELS 

Some  simple  AR  models  m=l  and  m=2  are  analyzed 
to  show  how  results  may  be  obtained.  Note  that  in 
all  fully  two-sided  AR  models  <p  ^ 

The  two-sided  simplest  AR  model,  m=l,  from 
(2.6)  is 

z  =  4>.  z  -+(}>  ,z  _+a  =<J>_  (a  ,+z  t,)+a  , 
x  T1  x— 1  -1  x+1  x  1  x-1  x+1  x 

al=(l-2o  4>  )-1  a2  .  (4.1) 

z  11  a 


^i-l^i+l1  for  a11  ^°*  then 

£  2 
P£apl'  <^i==p1/{1+p2)*Pl/a+Pl)  9  a  parabola.  (4.2) 

Note  1 l <% »  and  given  4>^,  p^=0.5 

2  L 

{^(1-4^)  }.  The  values  of  {^p^p^p.^}  are 
given  in  Table  2, 

The  AR  model  written  as  an  MA  model  is: 

v  {i!o  *i<vB;1,i}vv*i(wvi> 


+  + 


if  v  *-i  ax+i'  V1+*i(Vir-i) » 

V  VW'  W-l’+ll1,-24V  •••' 


*1  (1ri-l+1,l+l) '  <Tr-i-l+1t-i+l)  1 


and  tt  1*^1  • 

Given  <Plt  -  J0  <21)  4>2i  , 


_  ?  *2i-l 

”1  'i^l  (  i-1  »  *1  ' 


.  r  ,2i-l.  .21-1 

V  l  i-2  *1 

For  ^  =  .2 


TT  .  r  /  2i  .  .  2i 

W2  i=l  *i-l}  ^1  ' 


„  ®  ,  2i  %  .21 

VA  i-2  *i  •  “tc- 


IT  1  =  .22772201,  consequently  *  1.09108881, 


and  It^/ITq  =  .20871  *  ,  which  checks  Table  2. 

The  autoregressive  model  given  in  (4.1)  is 
simulated  given  that  0^  «  .2,  and  a^’s  *r*  N(0.1) . 

The  number  of  observations  in  each  simulation  is 
100.  First  a^’s  are  generated  and  z^'s  are 

obtained  by  using  (4.1):  to  do  this  consecutive 
forward  and  backward  substitutions  are  performed 
until  convergence  is  assured.  Results  of  the 
twenty-five  simulation  runs  are;  ■  .2282, 

=  0.0143,  B_  *  .2091,  of  ■  .0099.  For  one 
rl  1  91 

particular  run  when  0^  *  .2,  theoretically  P1  = 

.2087  and  0^  ^  n_1 (1-8?) /(l-p?)« (l-.04)/<100 <1- 
®1  1  1 

.0427))  =  .01.  Simulated  case  gives  r^  =  .1556, 
6^.1516,  Oq  %  n_1(l-e2)/(l-r^)  -  (1-.0231)/ 

(100 (1-. 0242) )  =  .01.  Since  using  $  =p 

a  w j  1  1^2^ 

(l+p0)  is  a  least  squares  estimate  so  072  = 

2  2  * 
tl-4^ ) /n (1-p^) ,  then  confidence  intervals  may  be 

found  assuming  the  a's  are  distributed  normally. 
The  minimum  variance  estimate  of  may  be 

obtained  as  indicated  in  the  case  of  0*s. 

Another  simple  two-sided  AR  model,  m=l,  is: 

zx  =  Wl+Z*+1>  +  *2(Zx-2+Zx.2>  +  *x  • 

°l  “  oju-ap^pj^"1  . 

P*  =  VWW*  *2(Pt-2+PU2)'t'10’  <4-4) 

the  Yule-Walker  equations  are: 

P1  *  <*>l(1+p2)  +  *2lPl+<V 

P2  “  *1(P1+P3J  +  <J>2(1+P4)  (4*5) 

♦j  =lPl<l+P4>  -  P2<P1+P3)^®1+P2)  U+P^-tPj+Pj)2] 
*2  -[P2(l+P2)-P1(P1+P3)V{(l+(»2)  (i+P4)-(P1+P3)2] 

(4.6) 

1 |  +  t  I <la '  a*14*  the  «<jviivalent  MA  model  is 

ifo  {^X+B;1>  +  \  •  <4. 7) 

Given  a  permissible  set  (P^Pj 'P3'P4} » 

is  determined.  Given  then  {p,  ,p„ ,p,  ,p^ } 

12  12  3  4 

is  determined  from  the  corresponding  MA  model, 
Aroian  and  Schmge  £1980) .  The  variances  and  the 
covariance  of  way  be  found  as  indicated  in 

Aroian  and  Schmee  (1980) ,  gnd^as  usual  the  mini¬ 
mum  variance  estimates  of  The  simula¬ 

tion  and  prediction  of  such  a  model  is  essential¬ 
ly  the  same  as  given  in  Box  and  Jenkins  (1976) 
for  pure  time  series,  but  four  starting  values 
will  be  needed  instead  of  only  two. 

For  m*2,  a  simple  two-sided  model  is: 

*x.,x  *  ^ll*x  -l,x  -l+^-l-llx,+l#*2+1 

1  (4.8) 

,x2-H+^-ll*x1+l  ,x2-l+*x1  ,x2 


with  solution: 


with  4>u  =  and  x 


*2  ' 


hcnco  z 


V 


-0.  (z 


1  X1_1,X2’1  X1+1,X2+1 
+<^2  <2xi-l  f  x2+l+ZXj  +  l  ,x2~l*  +a*i  '*2 

V(I-2puV2p.nV‘1  °a'  (4-9> 

pjj  j£  J  (p£  ^  ^  ^ 

4>2l0a  -l,*2+l+PJt  -1,1  +1)#  V  °'V°* 

The  characteristic  function  is: 


1-4  (B  B  ♦B'V1) 
1  x2  xA  x2 


f)  {B  B’1+B*1B  )  (4.10) 

1  X2  1  2 


which  may  be  used  to  find  the  corresponding  MA 
model  useful  to  obtain  all  needed  p?  p  ,  see 

1'  2 

Aroian  and  Schmee  (1980).  Note  |$.J  +  | $2 i <ls* 

The  Yule-WalJcer  equations  are: 

P11  ’  ^1(1+P22)  +  (^2(P02+P20) 

Pl-r  <*>l(P02+P20)  +  *2<1+P-22)'  ^  (4ai) 

(1+P-22)  ~P1-1  (P02+P2CP//^<1+P22)  '  (P02+P20)  ^ 

^2=*,t>l-l(1+P22)"Pll{P02+P2(j^^+P22)  ~(P02+P20)  ^ 

(4.12) 

which  involve  p1lfp.  and  p 


11 '  1-1  K02' K20 ' p2 , 2 


-2,2; 

A  a 


Given  estimates  of  these  six  correlations  ,4^) 
are  found  from  (4.12) . 

Conversely,  given  (<J> then  all  p0  0 

12  *i  #  X»2 

must  be  found  from  the  corresponding  MA  expansion. 
Estimation,  variances  and  covariance  of  ($^,$2) 

may  be  found  as  already  indicated  in  other  AR 
models,  particularly  the  minimum  variance  method. 
Another  obvious  two-sided  model  is: 

*  v  (*  +Z  )+<J>  U  T  X  +ZX  0.1  V  ) 

1#*2  A  X1,X2”"L  z  *2 

+  a  (4.13) 


V^OlV^loV'1  °a  ' 

Pi  H  £  +1^^2^PA  —1  &  +P£  +1  J,  ^ 

12  12  12  1212 

This  is  analyzed  in  exactly  the  same  way  as  the 
previous  models. 

Clearly,  I  ^  I  +  1  $2 1 <lj  $  the  characteristic 
function  is 


1  -  <f>  (B  +B'1)  -*2<B  +B"1)  , 

2  X2  1  X1 


(4.14) 


from  which  the  corresponding  MA  may  be  found  and 
all  pp  p  's.  The  Yule-Walker  equations  are: 

V*2 

P01  •  V1+P02>  + 

pio  *  Vpi-i+p-i-i)  *  V1+p20>'  <4-15) 


♦r  lp01  (1+p20>  -p10(pi-l+p-l-l)1/l  <1+P02»  <1+p20»  - 
(Pl-l+P-i -1,2]  ' 

*2tpio(1+po2,-piocpi-i+p-:.x'1/[(1+po2»<1+p2o>- 

(Pl-l+P-l-l)2j  ’  (4'X6) 

Estimation,  the  approximate  variance-covariance 
matrix,  and  confidence  limits  as  well  as  minimum 
variance  estimates  are  found  as  before. 

It  should  be  mentioned  that  the  partial  auto¬ 
correlation  function  of  the  AR  models  have  a  cut¬ 
off  property,  for  the  first  model  mwl,  all 
i>l  are  zero,  and  is  a  partial  coefficient  of 
correlation.  Where1there  are  two  $*s,  ^  and  <J>2 

both  are  partial  coefficients  of  correlation,  and 
4>i#  i  >  2  are  zero.  This  property  is  helpful  in 

determining  how  far  to  proceed  in  the  $*s;  alter¬ 
natives  are  the  analysis  of  variance  methods,  and 
that  in  which  the  variance  of  the  error  of  pred¬ 
iction  falls  as  the  #*s  increase.  Other  possible 
alternatives  still  unexplored  are  the  x2  test  and 
the  Kolmogoroff -Smirnov  test.  An  important  point 
to  remember  is,  as  the  number  of  $'s  increase,  so 
does  the  variance  of  the  forecast  errors.  * 


5.  ARMA  MODELS 

The  two  models  considered  are  simplifications 
of  (2.10)  and  C2.ll) : 


6,  (a 

1  x-1  x+1  x 


(5.1) 


*X1  ,X2  ,X2""^+Zxi  'X2*^  601Ux1,x2-i" 


a*l'X2+1*  +  aXl'X2 

(5.2) 

,1)  the  equations  for  o^,  p^. 

and  P2  are: 

a*  11+20^-40^) 

(5.3) 

i^U+P^V/ct*  {*107-261) 

(5.4) 

♦i  cpx+p3) + (6^-2*^)  o*/a* 

(5.5) 

From  these: 

al/al  “Ip2'^i(pl+p3>]/teJ-261i(i1]-Ip1-<|>1(l+P2)]/ 

[  {<(>10j-2e1>]-[i-2p1^/u+2e^-4e14llJ.  (s.o 

Thus  (5.6)  may  be  solved  for  using 

{p1#p^}.  The  restrictions  on  and 

^01  '^01^  *or  t*'e  ^  and  ^  apply.  Both  (5.1) 

and  (5.2)  may  be  written  in  the  equivalent  AR  or 
MA  model  since: 

{lHWC1,}‘*  ’  {1-VVBx1,}  v 

ax  -  <ifo  e}‘VBx1,1>tl-ei(VB»1}  *x  • 

and 

*x  “  {i£o  ^tVBxl)iHl_4,i<VBx1}  V  (5,7) 

with  similar  results  for  (5.2). 

A  more  general  model  than  (5.1)  is: 


j 


(5.8) 


z  «  <p  (z  +z  )  -  0  a  -0  a  +a  ,  (5.8) 

x  1  x— 1  x+1  1  x-1  -1  x+1  x 

which  reduces  to  (5.1)  if  0  ^  -  0^ . 

Now 

o^a+e2+02i.2^er2^e.1)  u^p^r1  , 

Pl=<ll  <1+p2* +  l°a/0z^  ^ieie-l”6l'0-l^  * 

P2=“>1  <P1+P3* +  [Ca/0z^  ^1  +ei6-l}  * 

p -<t>  (p-+p,)  +  [oWi  et>.e  e  )  .  (5.9) 

3124  a  Z  11-1 

2  2 

Solve  (5.9)  using  O^/a^  from  the  first  equation, 

substitute  this  into  the  other  three  and  solve 
for  using  the  sample  values 

{ r  ,r  ,r  ,r  }  for  the  p#s.  Write  (5.8)  as 
12  3  4 


{l-<fr.(B  +B*1)}  z  =  (1-6,  B  -0  _  B'1}  a  ,  (5. 

1  X  X  X  lx  —1  X  X 

then  z^  as  an  MA  process  is: 

z  =  (l-O.B  -6  , B-1} i l-<t>  (B  +B'1)}"1  a 

X  lx  —1  X  1  X  X  X 


10) 


z  =  [  { 1-0.  B  -0  b'TL  4l(B  +B  )x]  a  ,(5.11) 
x  lx  —1  x  i=0  1  x  x  x 

I  01 1  + 1  e_i | <x»  the  same  conditions  are  in 

(3.5)  and  (4.1) . 

Represent  as  an  AR  model  from  (5.10) 

a  =  {l-$_ (B  +B-1)}{l-e,B  -6  tB-1}"1  z 
x  Ixx  lx  -lx  x 


ax  Hi-^cvb;1,}  ii0{e1Bx+e_lB;1}i]zx  (5.12) 

with  the  same  restrictions  on  ,0  ,0  as  in 
(5.11). 

Now  suppose  0  ,6^,8  ^  are  given  satisfying  the 

restrictions  in  (5.11)  ,  what  are  the  values  of  p^, 
the  autocorrelation  function? 

Estimation  proceeds  as  indicated  in  Aroian  and 
Taneja  (1980) ,  by  changing  an  ARMA  model  to  an 
equivalent  AR  model  and  u9ing  the  results  from 
least  squares. 


6.  EXAMPLES 

Some  examples  will  be  completed,  particularly 
Whittle's  example  of  wheat  data  and  possibly 
some  others.  Simulations  will  be  done  in  a 
separate  study  as  well  as  further  extensions  of 
these  models.  One  sided  and  mixed  models  will  be 
done  in  the  future. 

7.  ISOTROPIC  PROCESSES 

Let  the  variable  x  represent  the  distance  from 
any  point  (X^X^  in  plane,  or  the  point 

( X i * X2 #X in  apace.  Then  for  any  stationary 

isotropic  process  the  results  from  Box  and 
Jenkins  (1976)  may  be  used  in  all  cases  for  MA, 
AR,  ARMA  replacing  t  by  x*  This  applies  not 
only  to  stationary  processes  but  to  nonstationary 
processes  if  one  uses  differences  as  indicated 
there.  Since  the  method  is  strightforward,  no 
further  discussion  is  needed.  For  isotropic 

processes  in  time  and  space,  in  place  of  z  .  as 

x,  t 

given  in  Aroian  et  al,  one  would  replace  x  by  X* 
and  retain  t  and  use  the  methods  indicated  there 


for  stationary  processes.  For  nonstationary 
processes  differences  in  two  variables  may  be 
used  or  transformations.  Another  alternative  to 
transformations  or  to  differencing  direct  treat¬ 
ment  of  nonstationarity  is  feasible  and  will  be 
investigated  subsequently. 

8.  CONCLUSIONS 

The  methods  of  time  series  in  m  dimensions  are 
applied  to  two-sided  spatial  models  in  one  and 
two  dimensions:  MA,  AR,  and  ARMA  models  illus¬ 
trate  the  techniques  including  estimation.  These 
results  presented  in  this  paper  are  based  on  the 
second  order  moments,  and  MA,  AR,  and  ARMA  models 
as  developed  in  time  series  in  m  dimensions. 

Papers  in  bibliography  numbered  as  1,2,3,4,11,12, 
14,  and  15  reflect  the  Aroian  point  of  view. 

Some  other  points  of  view  related  to  these  results 
may  be  found  in  Bartlett  (1975) ,  Bennett  (1979) , 
Besag  (1972) ,  Cliff  and  Ord  (1973)  and  Ord  (1975) . 
Bartlett  reflects  a  position  from  partial  differ¬ 
ential  equations,  and  power  spectrums  to  AR 
models,  a  broad  point  of  view  ewsring 
briefly  most  of  the  previous  work  before  1975. 
Bennett  covers  the  ideas  quite  thoroughly  and 
presents  a  comprehensive  bibliography,  but  does 
not  give  enough  details  as  Box  and  Jenkins  (1976) 
do  in  their  work.  Ord  considers  only  first  order 
autoregressive  models,  m=l ,  which  are  restricted 
and  not  general. 
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Table  2 


Values  of  (4^*  »  P2» 


*1 

1 

±  4 

P2 

P3 

0 

0 

0 

0 

±.05 

±.0501 

.0025 

±.0000 

±.l 

±.1010 

.0102 

±.0010 

±.2 

±.2087 

.0436 

±.0091 

±.25 

±.2679 

.0718 

±.0192 

±.3 

±.3333 

.1111 

±.0370 

±.4 

±.75 

.5625 

±.4291 

±.5 

±1 

1 
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